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Abstract 

We perform hydrodynamic derivations of the entrainment of matter (gas) induced by a 
surface elastic wave (representing possibly vacuum fluctuations) propagating along the flexi- 
ble vacuum-matter interface by considering the nonlinear coupling between the interface and 
the rarefaction effect via a boundary perturbation method. The critical reflux values asso- 
ciated with the product of the second-order (unit) body forcing and the Reynolds number 
(representing the viscous dissipations and is the ratio of wave inertia and viscous dissipation 
effects) decrease as the Knudsen number (representing the rarefaction measure) increases 
from zero to 0.1. We obtained the critical bounds for matter- freezed or zero- volume- flow- 
rate states corresponding to specific Reynolds numbers and wave numbers which might be 
linked to the dissipative evolution of the Universe. Our results also show that for certain 
time-averaged evolution of the matter (gas) there might be existence of negative-pressure 
states. 

Keywords : cosmic flows, physics of the early universe, gravity waves/theory 
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1 Introduction 

In today's standard inflationary scenario the amplified vacuum fluctuations of the inflation field 
provide the seeds for large-scale structure formation [1]. It is thus of interest to understand the 
effect of these quantum fluctuations on the universe dynamics and the implications they may 
have for cosmological observations. 

The mean cosmic density of dark matter (plus baryons) is now pinned down to be only ca. 
30% of the so-called critical density corresponding to a 'flat'-Universe. However, other re- 
cent evidence — microwave background anisotropics, complemented by data on distant super- 
novae — reveals that our Universe actually is 'flat', but that its dominant ingredient (ca. 70% 
of the total mass energy) is something quite unexpected: 'dark energy' pervading all space, 
with negative pressure. We do know that this material is very dark and that it dominates the 
internal kinematics, clustering properties and motions of galactic systems. Dark matter is com- 
monly associated to weakly interacting particles (WIMPs), and can be described as a fluid with 
vanishing pressure. It plays a crucial role in the formation and evolution of structure in the 
universe and it is unlikely that galaxies could have formed without its presence [1-2]. Analysis 
of cosmological mixed dark matter models in spatially flat Friedmann Universe with zero A term 
have been presented before. For example, we can start from the Einstein action describing the 
gravitational forces in the presence of the cosmological constant [3] 



where S rn is the contribution of the matter and radiation, kjj = 8^(7 at = 8ir/Mp, Gjy is the 
Newton's constant and Mp = 1.22 x 10 19 GeV/c 2 is the Planck mass. The set of equations 
governing the evolution of the universe is completed by the Friedman equations for the scale 
factor (a(t)) 
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where k = is for the flat-Universe, p and p is the density and pressure [4-5]. A large majority 
of dark energy models describes dark energy in terms of the equation of state (EOS) pd = to pj, 
where u is the parameter of the EOS, while pd and pd denote the pressure and the energy density 
of dark energy, respectively. The value to = — 1 is characteristic of the cosmological constant, 
while the dynamical models of dark energy generally have to > —1. 

Meanwhile, it is convenient to express the mean densities pi of various quantities in the Universe 
in terms of their fractions relative to the critical density: Qi = Pi/p cr it- The theory of cosmo- 
logical inflation strongly suggests that the total density should be very close to the critical one 
{Vttot ~ 1), and this is supported by the available data on the cosmic microwave background 
(CMB) radiation. The fluctuations observed in the CMB at a level ca. 1CT 5 in amplitude ex- 
hibit a peak at a partial wave I ~ 200, as would be produced by acoustic oscillations in a flat 
Universe with Q, to t ~ 1- At lower partial waves, / » 200, the CMB fluctuations are believed to 
be dominated by the Sachs- Wolfe effect due to the gravitational potential, and more acoustic 
oscillations are expected at I > 200, whose relative heights depend on the baryon density f2&. At 
even larger values, / > 1000, these oscillations should be progressively damped away [1-2,3,6]. 
Influential only over the largest of scales-the cosmological horizon-is the outermost species of 
invisible matter: the vacuum energy (also known by such names as dark energy, quintessence, x- 
matter, the zero-point field, and the cosmological constant A). If there is no exchange of energy 
between vacuum and matter components, the requirement of general covariance implies the time 
dependence of the gravitational constant G. Thus, it is interesting to look at the interacting 
behavior between the vacuum (energy) and the matter from the macroscopic point of view. One 
related issue, say, is about the dissipative matter of the flat Universe immersed in vacua [7] and 
the other one is the macroscopic Casimir effect with the deformed boundaries [8]. 
Theoretical (using the Boltzmann equation) and experimental studies of interphase nonlocal 
transport phenomena which appear as a result of a different type of nonequilibrium represent- 
ing propagation of a surface elastic wave have been performed since late 1980s [9-10]. These 
are relevant to rarefied gases (RG) flowing along deformable elastic slabs with the dominated 
parameter being the Knudsen number (Kn = mean-free-path/L^, mean-free-path (mfp) is the 
mean free path of the gas, Ld is proportional to the distance between two slabs) [11-13]. The role 
of the Knudsen number is similar to that of the Navier slip parameter N s [14]; here, N s = pS/d 
is the dimensionless Navier slip parameter; S is a proportionality constant as u s = St, t : the 
shear stress of the bulk velocity; u s : the dimensional slip velocity; for a no-slip case, S = 0, but 
for a no-stress condition. S = oo, p is the fluid viscosity, d is one half of the distance between 
upper and lower slabs). 

1.1 Present Approaches and Objectives 

Here, the transport driven by the wavy elastic vacuum-matter interface will be presented. The 
flat-Universe is presumed and the corresponding matter is immersed in vacua with the inter- 
face being flat-plane like. We adopt the macroscopic or hydrodynamical approach and simplify 
the original system of equations (related to the momentum and mass transport) to one single 



4 



higher-order quasi-linear partial differential equation in terms of the unknown stream function. 
In this study, we shall assume that the Mach number Ma <C 1, and the governing equations are 
the incompressible Navier-Stokes equations which are associated with the relaxed slip velocity 
boundary conditions along the interfaces [11-14]. We then introduce the perturbation tech- 
nique so that we can solve the related boundary value problem approximately. To consider the 
originally quiescent gas for simplicity, due to the difficulty in solving a fourth-order quasi-linear 
complex ordinary differential equation (when the wavy boundary condition are imposed), we can 
finally get an analytically perturbed solution and calculate those physical quantities we have in- 
terests, like, time-averaged transport or entrainment, perturbed velocity functions, critical unit 
body forcing corresponding to the freezed or zero-volume-flow-rate states. These results might 
be closely linked to the vacuum-matter interactions (say, macroscopic Casimir effects) and the 
evolution of the Universe (as mentioned above : the critical density [1-2]). Our results also show 
that for certain time-averaged evolution of the matter (the maximum speed of the matter (gas) 
appears at the center-line) there might be existence of negative-pressure states. 

2 Formulations 

2.1 Preview of Mesoscopic Approach 

The effects of elastic or deformable interfaces, like surface elastic waves (SLW) interacted with 
volume and surface phonons (propagating along the elastic boundaries), upon the entrained 
transport of gases have been studied since early 1950s [9-10], however, the mathematical difficulty 
is essential therein. The role of elastic macroscopic walls resembles that of microscopic phonons. 
As presented in [9], for the description of the transport processes in the non-equilibrium gas-solid 
system including the processes occurring in the case of propagation of an elastic wave in a solid, 
we need to solve 

g/(x,v,t) | df(x,v,t) _ 

dt + dx -wn), 

dn(x,kj,t) | dn(x,kj,t) _ 

dt + r ax " 4iWJ ' 

and the associated boundary conditions (scattering and interacting laws near the interface) 

\v r \f + = / dVi|«i|/~(Vi)W(Vi -► v r ), 
Jvi<0 

^n-(k,) = ]T [^n + (k n ) + N g (k n ) + N p (k n )]V p (k n - kj;co). 

1 k,i(ci>0) * 

/, n, and H denote the distribution function for gases, volume phonons, and surface phonons, 
respectively. I g , I v , and I s are the corresponding collision integrals. Please see the details in [9] 
for other notations or symbols. 
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As the details of the vacuum energy (also known by such names as dark energy, quintessence, 
the zero-point field, and the cosmological constant [4]) and dark matter (possible candidates : 
non-luminous baryonic matter in the form of large planets or dead stars or weakly interacting 
elementary particles that pervade large regions of space [2]) are not known so far and we only 
consider the overall effect of the matter, thus the detailed microscopic collision or scattering oc- 
curring along the interface and transfer of the energy from matter to the vacuum energy (which 
could be cosmological constant) will not be treated here and is beyond our present approach. 
We noticed that particle physics candidates for dark matter fall into three basic categories de- 
pending upon their masses and interactions [2]. Weakly interacting particles that are in thermal 
equilibrium at a very early stage in the evolution of the universe will eventually fall out of 
equilibrium as the universe expands. The decoupling time (or temperature) when this occurs 
depends on the expansion rate of the universe as well as the couplings of these particles to 
other particles that are still in equilibrium. Particles that are (non-)relativistic at the time that 
galaxies start to form are referred to as (cold) hot dark matter. The simplest example of hot 
dark matter is a neutrino with a very small mass (< O(20) eV) and of cold dark matter a very 
heavy neutrino (~ 100 GeV). In both cases the interaction rates are determined by the standard 
model of electroweak interactions (SM). The third type of particle dark matter can arise during 
the QCD phase transition as the universe cools down. In this case the result can be a gas of 
axions [4].) 

To escape from above (many-body problems) difficulties, we plan to use the macroscopic ap- 
proach which is a complicated extension of previous approaches [10] because we still take the 
rarefaction effect (due to the incomplete accomodations of the momenta or energy for colli- 
sions or reflections of particles from the matter-vacuum interface which are approximated and 
represented by the Knudsen number or the slip velocities [11-13]) into account. 

2.2 Interface Treatment 

We consider a two-dimensional matter-region of uniform thickness which is approximated by a 
homogeneous rarefied gas (Newtonian viscous fluid). The equation of motion is 

d 2 u 

(Ai + ^)graddivu + ^Vu + pp = /9_-, (1) 

where and fi are Lame constants, u is the displacement field (vector), p is the mass den- 
sity and p is the body force for unit mass. The Navier-Stokes equations, valid for Newtonian 
fluids (both gases and liquids), has been a mixture of continuum fluid mechanics ever since 
1845 following the seemingly definitive work of Stokes and others [15], who proposed the fol- 
lowing rheological constitutive expression for the fluid deviatoric or viscous stress (tensor) T : 
T = 2//Vu + X L IV ■ u. 

The flat-plane boundaries of this matter-region or the vacuum-matter interfaces are rather flex- 
ible and presumed to be elastic, on which are imposed traveling sinusoidal waves of small 
amplitude a (possibly due to vacuum fluctuations). The vertical displacements of the upper 
and lower interfaces (y = d and —d) are thus presumed to be rj and — rj, respectively, where 
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rj = acos[2-7r(x — ct)/X], A is the wave length, and c the wave speed, x and y are Cartesian 
coordinates, with x measured in the direction of wave propagation and y measured in the direc- 
tion normal to the mean position of the vacuum-matter interfaces. The schematic plot of above 
features is shown in Fig. 1. 



a = ed 


y 

r/(x, t) = a cos ^(x — ct) 
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wave speed c 



Fig. 1 Schematic diagram of the wavy motion of the vacuum-matter interfaces. 

It would be expedient to simplify these equations by introducing dimensionless variables. We 
have a characteristic velocity c and three characteristic lengths a, A, and h. The following 
variables based on c and d could thus be introduced : 



, x , y , u , v , rj , ip , ct , p 

v = 3> y = 3' u = -■> v = V = j, w = — -,, t = —, p = 

a a c c a cd d pc z 



where ip is the dimensional stream function, u and v are the velocities along the x- and y- 
directions; p is the density, p (its gradient) is related to the (unit) body forcing. The primes 
could be dropped in the following. The amplitude ratio e, the wave number a, and the Reynolds 
number (ratio of wave inertia and viscous dissipation effects) Re are defined by 

a 2nd cd 

e=-, a = — — , Re= — . 
a A v 

We shall seek a solution in the form of a series in the parameter e : 

ip = tpo + eVd + e 2 V>2 H , 

dp dp dp 2 dp 

with u = dtp/dy, v = —dip/dx. The 2D (x- and y-) momentum equations and the equation 
of continuity could be in terms of the stream function ip if the p-term (the specific body force 
density, assumed to be conservative and hence expressed as the gradient of a time-independent 
potential energy function) is eliminated. The final governing equation is 

|v^ + ^fc-feV^ = ^VV, + (2) 

and subscripts indicate the partial differentiation. Thus, we have 
d 1 

— V 2 V>o + V'OyVVoz - VfeVVch/ = ^ v4 ^o, (3) 



— V 2 ^l + TpOyV 2 7pl x + lplyV 2 7p 0x ~ VfeV^ly - lpl x V 2 7pOy = ( 4 ) 

V'OxV 2 ^ - Vi^V 2 ^^ - feV 2 ^ = ^ v4 ^2, (5) 

and other higher order terms. The (matter) gas is subjected to boundary conditions imposed by 
the symmetric motion of the vacuum-matter interfaces and the non-zero slip velocity : u = =F 
Kn du/dy [11-13], v = ±drj/dt at y = ±(1 + rj), here Kn=mfp/(2d). The boundary conditions 
may be expanded in powers of r/ and then e : 

tpoyh + e[cosa(x - t)ip 0yy \i + V'lj/li] + e 2 ^ ^ 1 cos 2 a(x - t) + ^2y\i+ 

cosq(x - t)ip lyy \i] H = -Kn{^ ?/y|i + e[cosa(x - i)Voj/ W li + 

e 2 ^0yy W |l cQg 2 ^ _ ^ + COS a(x - + V^j/ll] H }, (6) 

^Ox|i + e[cosa(x - t)ifj 0xy \i + ip lx \i] + e 2 [^2^!i C os 2 a(x - t)+ 
cosa(x — t)tpi xy \i + ^>2x|i] + • • • = — easina(x — t). (7) 
Equations above, together with the condition of symmetry and a uniform (dp/dx)o, yield : 

^ = tf [(l + 2Kn)y-^], K = ^("^)o, (8) 

V'i = \{4>{yy a{x ~ t] + <f (y)^^}, (9) 

where the asterisk denotes the complex conjugate. A substitution of ipi into Eqn. (3) yields 

d 2 d 2 
{— j - a 2 + i<xRe[l - if (l - y 2 + 2Kn)]}(— , - a 2 )<j) - 2iaK Re(j) = 
dy z ' dy z 

or if originally the (matter) gas is quiescent : K = (this corresponds to a free (vacuum) 
pumping case) 

d 2 d 2 

(to -" 2 )(-r^ -a 2 U = 0, a 2 = a 2 -i<xRe. (10) 

ay" 1 ay z 

The boundary conditions are 

^(±l)±^ w (±l)Kn = 2^ (l±Kn) = 0, 0(±1) = ±1. (11) 

Similarly, with 

y, 2 = i{D(y) + E{y)e l2a ^ + ^(y)^ 2 ^-*)}, (12) 



we have 



D yyyy = -^o~(H>ly ~ 4>*<l>yy)y, (13) 
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d 2 d 2 
[—= - (4a 2 - 2iaRe)}(—^ - Aa 2 )E - i2aReK (l - y 2 + 2Kn) 
dy z dy z 

dp* zotRs 
(^2 " 4a ^ E + ^K ReE + _^_(^ w - = 0; (14) 

and the boundary conditions 

D y (±l) + \[<P vy (±l) + ^(±1)] - 2K = TKn{i[^ ra (±l) + <f>* yyy (±l)} + D n (±l)}, (15) 

E V {±1) + ^ w (±l) - ^ = TKn[^ yra (±l) + E yy (±l)], (16) 

£(±1) + ^„(±1) = (17) 
where ifo is zero in Eqns. (13-16). 
2.3 Originally Quiescent Gas 

To simplify the approach and obtain preliminary analytical solutions of above complicated equa- 
tions and boundary conditions, we only consider the case in which (dp/dx)o or Kq = ipQ = 0. 
After lengthy algebraic manipulations, we obtain 

cj) = CQ e ay + cie~ ay + c 2 e ay + c 3 e"^, 

where c = (A + A )/Det, c x = -{B + Bo)/ Deb, c 2 = (C + C )/Det, c 3 = — (T + T )/Det; 

Det = Ae a - Be~ a + Ce a - Te~ a , 

,4 = e Q a 2 (r 2 e- 2a - s 2 e 2a ) - 2aae- Q u; + aae Q z( e - 2a r + e m s), 
A = e- Q a 2 (r 2 e - 2a - s 2 e 2a ) + 2aae Q z - aae" a u;(e 2a s + e" 2a r), 
5 = e - Q Q 2 (r 2 e- 2a - s 2 e 2a ) + 2aae a z - aae- a w{e~ m r + e 2a s), 
B = e Q a 2 (r 2 e- 2a - s 2 e 2a ) - 2aae~ a w + aae a z{e~ m r + e 2a s), 
C = e~ a aa(wse a ~ a — rze a ~ a ) — ae 2a+a z(az — as) + ae~ a w(ae a ~ a w — ae a ~ a r), 
C = e a aa{wse«- a - rze a ~«) - az(zae 2a ~« - ae a s) + aw{ae-^ +2a) w - ae- (2o+a) r), 
T = e- a aa(zse a+a - ru> e - (a+a) ) - aa(e 2a " a zr - e a ws) + a 2 e - a (-e 2a z 2 + e~ 2a w 2 ), 
T = e a aa(zse" +a - rwe~ {a+5!) ) - aa(e~ a zr - e*~ 2a ws) + a 2 e a {-e 2a z 2 + e~ 2a w 2 ), 
with r = (1 — «Kn), s = (1 + «Kn), w = (1 — <xKn), z = (1 + aKn). 

To obtain a simple solution which relates to the mean transport so long as only terms of 0(e 2 ) 
are concerned, we see that if every term in the x-momentum equation is averaged over an interval 
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of time equal to the period of oscillation (due to vacuum fluctuations) , we obtain for our solution 
as given by above equations the time-averaged (unit) body forcing 



where IIo is the integration constant for the integration of equation (12) and could be fixed 
indirectly in the coming equation (22). Now, from Eqn. (14), we have 

£> y (±l)±Kn£> w (±l) = -^[<M±1) + ^(±1)] TKn{l[<^(±l) + <^(±1)]}, (19) 
where D y (y) = Hoy 2 + a\y + 02 + C(y), and together from equation (12), we obtain 

C(y) = \ C ^ C (^*)V + C *0 C 2 c (a+a)y + C C 3 c (a-a*)y + C C 3 e (q-5)y , 

2 g\ 92 $ 9l 

^C2 e (a*~a)y + ^_ e (.a-a)y + £l£3 e -(a*+a)j/ + ^£3_ e -(a+a)y + 
ffl 54 5l 52 

^3 e (a-a*)y + 4^3 e (<5*-a)j/ + 2 ^c| e (a*+a)j/ + o C3 °3 c -(a*+a) M l ^0) 
ffi 5| 56 56 

with gi = a + a* , g2 = a + a, g% = a — a* , = a — a, g§ = a — a* , gQ = a + a* . In 
realistic applications we must determine IIo from considerations of conditions at the ends of the 
matter-region. a\ equals to zero because of the symmetry of boundary conditions. 
Once IIo is specified, our solution for the mean speed (u averaged over time) of matter-flow is 

U = e 2 ^ = j{C(y) - C(l) + R - KnC„(l) + U [y 2 - (1 + 2Kn)]} (21) 

where Ro = —{[4>yyiX) + $^(1)] — Kn[0 3/w (l) +(j)* yy (l)]} /2, which has a numerical value about 
3 for a wide range of a and Re (playing the role of viscous dissipations) when Kn= 0. To 
illustrate our results clearly, we adopt U (Y) = u(y) for the time-averaged results with y = Y in 
the following. 

3 Results and Discussion 

We check our approach firstly by examining Ro with that of no-slip (Kn= 0) approach. This 
can be done easily once we consider terms of D y (y) and C(y) because to evaluate Rq we shall at 
most take into account the higher derivatives of <f>(y), like <j> yy (y), 4>yyy{y) instead of <fi y (y) and 
escape from the prescribing of 02- 

Our numerical calculations confirm that the mean streamwise velocity distribution (averaged 
over time) due to the induced motion by the wavy elastic vacuum-matter interface in the case of 
free (vacuum) pumping is dominated by Ro (or Kn) and the parabolic distribution — Ho(l — y 2 ). 
Rq which defines the boundary value of D y has its origin in the y-gradient of the first-order 
streamwise velocity distribution, as can be seen in Eqn. (14). 

In addition to the terms mentioned above, there is a perturbation term which varies across the 
channel : C(y) — C(l). Let us define it to be 

F(y) = a^ [C(y) - C(1)] (22) 
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To compare with no-slip (Kn= 0) results, we plot three cases, a = 0.1, 0.4, and 0.8 for the same 
Reynolds number Re = 1 of our results : Kn= 0.1 with those Kn= into Fig. 2. We remind 
the readers that the Reynolds number here is based on the wave speed. This figure confirms 
our approaches since we can recover no-slip results by checking curves of Kn= and finding 
them being almost completely matched in-between. The physical trend herein is also the same 
as those reported in Refs. [12-13] for the slip-flow effects. The slip produces decoupling with 
the inertia of the wavy interface. 

Now, let us define a critical reflux condition as one for which the mean velocity U(Y) equals to 
zero at the center-line Y = (cf. Fig. 2). With equations (12,20-21), we have 



dp [ a 2 fe 2 F(0)/20Q + KnC'(l)-flo] 
n 0cr = «e(-) 2 = _ (1 + 2Kn ) (23) 

which means the critical reflux condition is reached when Ho has above value. Pumping against a 
positive (unit) body forcing greater than the critical value would result in a backward transport 
(reflux) in the central region of the stream. This critical value depends on a, Re, and Kn. There 
will be no reflux if the (unit) body forcing or pressure gradient is smaller than this IIo. Thus, 
for some IIo values less than IIo cr , the matter (flow) will keep moving or evolving forward. On 
the contrary, parts of the matter (flow) will move or evolve backward if IIo > no cr . This result 
could be similar to that in [16] using different approach or qualitatively related to that of [5] : 
even for very slow growth of A, the gravitationally bound systems become unbound while the 
nongravitationally bound systems remain bound for certain parameters defined in [5] (e.g., rj). 

3.1 Possible Massless and Negative-Pressure States 

We present some of the values of IIo (a, -^ e i Kn = 0, 0.1) corresponding ot freezed or zero- volume- 
flow-rate states (j\ U (Y)dY = 0) in Table 1 where the wave number (a) has the range between 
0.20 and 0.80; the Reynolds number (Re)= 0.1, 1, 10, 100. We observe that as Kn increases from 
zero to 0.1, the critical IIo ° r time-averaged (unit) body forcing decreases significantly. For 
the same Kn, once Re is larger than 10, critical reflux values IIo drop rapidly and the wave- 
modulation effect (due to a) appears. The latter observation might be interpreted as the strong 
coupling between the vacuum-matter interface and the inertia of the streaming matter-flow. The 
illustration of the velocity fields for those zero-flux (zero-volume-flow-rate) or freezed states or 
massless states are shown in Figure 3. There are three wave numbers : a = 0.2, 0.5, 0.8. The 
Reynolds number is 10. Both no-slip and slip (Kn=0.1) cases are presented. The arrows for slip 
cases are schematic and represent the direction of positive and negative velocity fields. 
Some remarks could be made about these states : the matter or universe being freezed in the 
time-averaged sense for specific dissipations (in terms of Reynolds number which is the ratio 
of wave-inertia and viscous dissipation effects) and wave numbers (due to the wavy vacuum- 
interface or vacuum fluctuations) for either no-slip and slip cases. This particular result might 
also be related to a changing cosmological term (growing or decaying slowly) or the critical 
density mentioned in Refs. [1-2]. If we treat the (unit) body forcing as the pressure gradient, 
then for the same transport direction (say, positive x-direction), the negative pressure (either 
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downsdtream or upstream) will, at least, occur once the time-averaged flow (the maximum 
speed of the matter (gas) appears at the center-line) is moving forward! For example, for 
Ho = Re{^)2 = —10 < (Re = l,a = 0.5), the velocity field (profile) is shown in Fig. 4. One 
possible p-pair for uniform (negative) gradient ( mean value theorem) 1 Pdownstream — Pupstream ^ 0, 
with X downstream Xupstream ^ . Pdownstream ^ 0, Pdownstream *^ Pupstream *^ 0. 

The possible explanation for this latter result (negative-pressure states) is : Please be reminded 
that our results are time-averaged (averaged over an interval of time equal to the period of 
oscillation (due to vacuum fluctuations)) and are already in steady state or the 2nd.-order 
transport (even though the primary state of the matter gas is quiescent without the Big-Bang 
explosion-driven force) is almost fully developed (which means the viscous dissipation being 
balanced by the wave inertia coming from the elastic interface or vacuum fluctuations). The 
quasi-stationary state after taking time-average is a balanced combination of both energy-input 
and energy-dissipation with the former coming from the small-amplitude elastic deformations of 
the entire configuration (elastic energy) and the latter coming from the incomplete exchange of 
the momenta or energy along the vacuum-matter interface and the internal friction of the matter 
(Navier-Stokes gases with viscosities prescribed). The present approach is not exactly the same 
as previous approaches (say, [5,9]). There are indeed perturbed 2nd. -order time-averaged effects 
and parts of them could induce the negative-pressure states demonstared above (cf. Fig. 4) for 
specific Reynolds numbers, wave numbers, and Knudsen numbers (one is illustrated above)! 
Under the same situation, we further demonstrate the effects of viscous dissipation in Figs. 5 
(Re = 100, a = 0.2, n = -50) and 6 (Re = 0.1, a = 0.2, n = -0.1). Once d and c (the 
wave speed) are fixed, the former corresponds to the case of less viscous dissipation while the 
latter corresponds to the case of much more viscous dissioation. Meanwhile, the time-averaged 
transport induced by the wavy interface is proportional to the square of the amplitude ratio 
(although the small amplitude waves being presumed), as can be seen in Eqn. (12) or (20), 
which is qualitatively the same as that presented in [9] for analogous interfacial problems. In 
brief summary, using the boundary perturbation methods, the entrained transport or pattern 
: either postive or negative and there is possibility : freezing (which could be linked to the 
dissipative evolution of the Universe) due to the wavy vacuum-matter interface is mainly tuned 
by the (unit) body forcing or Ho for fixed Re. Meanwhile, Yi.Q cr depends strongly on the Knudsen 
number (Kn, a rarefaction measure) instead of Re or a. We also noticed that it has been recently 
pointed out that fluctuations in the stress-energy tensor of quantum fields may be important for 
some states in curved spacetimes or even flat spacetimes with nontrivial topology [16]. Hu and 
Phillips, for instance, computed the energy density fluctuations of quantum states in spatially 
closed Friedmann-Robertson- Walker models and showed that these could be as important as the 
energy density itself [17]. If so, these fluctuations may induce relevant backreaction effects on the 
gravitational field (the spacetime geometry). Some of our results shown above resemble these 
latter results. We hope that in the future we can investigate other issues like the role of phase 
transition and that of cyclic universes [18-21] using the present or more advanced approach. 
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Fig. 2 Demonstration of Kn and a effects on the time-averaged velocity (fields) profiles. 
ITo = IIo cr . The Reynolds number (the ratio of the wave- inertia and viscous dissipation 
effects) is 1. U(Y) = at Y = 0. Kn is the Knudsen number (a rarefaction measure). 
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Fig. 3 Demonstration of the zero-flux states : the mean velocity field U (Y) for 

wave numbers a = 0.2, 0.5, 0.8. The Reynolds number is 10. Kn is the rarefaction measure 

(the mean free path of the particles divided by the characteristic length) . 

The arrows are schematic and illustrate the directions of positive and negative U(Y). 

The integration of U(Y) w.r.t. Y for these velocity fields gives zero volume flow rate. 
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Fig. 4 Demonstration of the negative-p states : the mean velocity field U (Y) 
for the wave number a = 0.2 and the Reynolds number Re = 10. 
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Fig. 5 Demonstration of the negative-p states : the mean velocity field U (Y) 
for the wave number a = 0.2 and the Reynolds number Re = 100. IIo = —50. 
As Re = c d/u, once c (the wave speed) is fixed, this transport corresponds to 
the case with less viscous dissipation. 



18 



=) 



1.7 
1.6 
1.5' 
1.4 
1.3 
1.2 
1.1 



Re=0.10, n Q =-0.10, a=0.20 



~s~ Kn=0 
-x- Kn=0.1 



-0.5 




Y 



0.5 



Fig. 6 Demonstration of the negative-p states : the mean velocity field U (Y) 
for the wave number a = 0.2 and the Reynolds number i?e = 0.1. IIo = — 0.1. 
As Re = c d/v, once c (the wave speed) is fixed, this transport corresponds to 
the case with much more viscous dissipation. 



Table 1: Zero-flux or freezed states values (IIq) for a flat vacuum-matter interface. 







Re 








Kn 


a 


0.1 


1 


10 


100 





0.2 


4.5269 


4.5269 


4.5231 


4.3275 


0.5 


4.6586 


4.6584 


4.6359 


4.2682 


0.8 


4.9238 


4.9234 


4.8708 


4.4488 


0.1 


0.2 


2.4003 


2.4000 


2.3774 


1.2217 


0.5 


2.4149 


2.4132 


2.2731 


-0.9054 


0.8 


2.4422 


2.4379 


2.0718 


-3.4151 



